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Abstract 

A gauge- field theory for massive neutral particles is developed on the basis of the 
'53 ■ real four-component Majorana equation. By use of its spin operator, a purely 

^~>| imaginary representation of the SU(2) algebra can be defined, which gives a 

covariant derivative that is real. Such a coupling to the gauge field preserves 
the real nature of the Majorana equation even when including interactions. 
As the associated isospin is four-dimensional, this procedure introduces four 
intrinsic degrees of freedom to the Majorana field, which may be related to 
four flavours. The main aim is to describe here the mathematical possibility for 
coupling Majorana particles with a gauge field which resembles that of the weak 
interaction. By adding a fourth member to the family, flavour could become a 
dynamic trait of the neutral Majorana particles, and thus lead to a dynamic 
understanding of mixing. 
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1 Introduction 



According to the canonical standard model of elementary particle physics, the 
leptons and quarks come in three flavours, are massless and thus obey chiral 
symmetry, but then they can acquire mass through the Higgs mechanism (see, 
e.g., the text books [TJ [5J |3]). The Dirac equation [?] is fundamental in all 
this and well understood, however the nature of the neutrinos involved remains 
less clear. Are they Dirac fermions or Majorana [5] particles? The physics of 
neutrinos and its application to particle astrophysics remains a very active re- 
search area. In the past neutrinos were often described by the massless Weyl [6] 
equations involving only two-component Pauli [7] spinors. However, since con- 
vincing empirical evidence [8] for the finite neutrino masses and their associated 
oscillations [HI [TU] had been found in the past decades, massive neutrinos have 
been discussed, and furthermore another very massive neutrino species or sterile 
ones have been considered in four-neutrino models for example to explain 
the light neutrinos masses by the see-saw mechanism [9] . The tensions with the 
three- neutrino paradigm have recently been discussed by Kayser |12j , who also 
reviewed the key arguments in support for neutrinos being Majorana fermions 
[13] . Therefore the Majorana equation with various mass terms has found strong 
attention, either in its complex two-component (see the recent review by Dreiner 
et al. [14] ) or real four-component form, and been used in modern quantum field 
theory for the description of massive neutrinos. Thus theoretical reasoning and 
new empirical results from laboratory as well as cosmology suggest the possible 
existence of a fourth neutrino flavour. The state of affairs (as of 2006) and the 
research perspectives are described in the comprehensive review by Mohapatra 
and Smirnov [15] . 

The purpose of the present paper is to show that for massive neutral Majo- 
rana particles (perhaps representing the observed neutrinos) a gauge-field theory 
based on the real four-component Majorana equation is feasible. Using its spin 
operator, which obeys the angular momentum algebra like SU(2) but is purely 
imaginary, we can define an appropriate symmetry group which gives a connec- 
tion that is real, and thus the coupling to the gauge field (also real) does keep 
the real nature of the Majorana equation even when including interactions. As 
the isospin is four-dimensional, this procedure introduces four intrinsic degrees 
of freedom to the Majorana field, the interpretation of which remains open to 
speculation. Here the main aim is to describe the mathematical possibility for 
coupling the Majorana particles with a gauge field, which resembles that of the 
weak interactions between leptons and quarks [Jj)J H3 HE] ■ Although the neu- 
trinos are presently known to come in three flavours, we will here argue that, 
by adding a fourth member to the family, flavour could become a dynamic trait 
of the neutral Majorana particles and be related to the gauge theory described 
subsequently. 

The paper is organized as follows. We discuss the relevant key aspects of the 
Majorana equation (its eigenfunctions are given in the appendix) in the next 
section, which is more of a tutorial nature. Then the real Lorentz transfor- 
mation is given, in which the purely imaginary Majorana spin operator occurs 
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prominently. In our opinion, it strongly suggests itself as the adequate isospin 
for the SU(2)-like algebra adopted in the gauge symmetry considered in the 
subsequent section. When assuming four flavours, the resulting gauge theory 
provides a dynamic picture of flavour mixing, i.e. quantum flavour dynamics, 
and likely also new insights into neutrino oscillations. 



2 The Majorana equation 

In this introductory section, we consider the Dirac equation in its Majorana 
representation. We use standard symbols, notations and definitions, and con- 
ventionally units of h = c = 1, with the four- momentum operator denoted as 
= (E, — p) = idp = i(d/dt,d/dx), which acts on the spinor wave function 
ip(x,t). The particle mass is m. The four- vector 7^ consists of the four Dirac 
gamma matrices that come in various representations. As first found out by Ma- 
jorana [5 , there exists a purely imaginary representation such that 7^ = ry^ 1 , 
which makes the Dirac equation real. We simply refer to it as Majorana equation 
in the remainder of this work. This real equation reads 

7"o^ + rmp = 0. (1) 

The Majorana equation involves the below defined 4x4 matrices. Through- 
out the text we will use top-barred symbols to indicate real 4x4 matrices. The 
gamma matrices can after [19 be defined as 

The gamma matrices mutually anticommute and obey: 7^7^ + 7^7^ = —2g^, 
with euclidian metric g^ v . The three associated real 2x2 matrices read 

1(M „ ( 1 A ( -1 

-1 



1 M= 1 • ( 3 ) 



Note that 7 = /3a, and that all three matrices mutually anti-commute with 
each other, just like the Pauli matrices do. All what we need in the following 
are the algebraic properties of these 2x2 matrices involved, like a/3 + /3a = 0, 
«7 + 7a = 0, and j3j + 7/3 = 0, and that a 2 — (3 2 — 1 and 7 2 = — 1. The three 
Pauli matrices have their standard form, and are given by cr x = j3, a y — i7,and 
er z = a. 

The four-component Majorana equation ([l} is a direct consequence of, and 
fully equivalent to, the two-component complex Majorana equation. It involves 
only the Pauli matrix operators acting on a two-component spinor (f>, but in- 
troduces subtle complications that are caused by the spin-flip operator r, and 
reads as follows: 
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It can be derived without invoking the Dirac equation in the first place, as was 
demonstrated by Marsch [19l [20] recently. As shown years ago by Case [21] , 
this equation can be also derived from Dirac's equation in its chiral form. 

Above in (Q} we defined an important operator that is not of pure alge- 
braic nature but involves the complex-conjugation operator named as C. This 
antihermitian operator called r is defined as r = a y C, and obeys r _1 = — r. 
The operation of r on the spin vector er leads to its inversion, i.e., the oper- 
ation Terr -1 = <x yields a spin flip. We also have ri = — ir because of the 
action of C. Therefore, r anti-commutes with the momentum four-vector oper- 
ator P M . Moreover, a simple phase factor like exp(i0) (with some angle 9) does 
not commute with the mass term on the right-hand side of (U]). Consequently, 
that equation can neither describe electromagnetic interactions nor most other 
complex gauge-field couplings in which the imaginary unit i appears explicitly. 

To avoid the complications introduced by the operator r, we prefer to work 
with the real four-dimensional Majorana equation ([!}, which we quote again 
but now in conventional Hamiltonian form as used by Majorana [5] himself as 
follows 

d d 



ill ' a ' dx) ^ x '*) = n^fat), (5) 

which formally resembles the complex two-component version (U). This last 
equation is as usually obtained by multiplying the manifestly covariant form |T]) 
from the left by 7 = /? and defining <5 M = — = (l,ck). The real matrix 
three-vector a is then given as 



P \ f 7 \ f a 
ft J ' 1 -7 J ' 1 a 



(G) 



and is symmetric, which means equal to its transposed matrix, a = ct T . The 
eigenfunctions of the Majorana equations ([4]) and (0 are briefly discussed in 
the appendix Section [6] where it is shown that the free real Majorana can be 
decomposed into left- and right- helical particle and antiparticle components. It 
is chirally irreducible, though, unless one considers complex solutions. 

At this point, we quote the Lagrangian density of the Majorana field, which 
is obtained by inserting the Majorana gamma matrices into the general Dirac 
Lagrange density. Note that the factor i is important here to ensure that we 
are dealing with hermitian operators having real eigenvalues. Yet the factor i is 
irrelevant, and does not appear, in the Majorara equation ([5]) itself. The final 
result is (where the superscript T denotes the transposed spinor) given by 

C M =i^ T (a M ^-^m)^. (7) 

Its variation with respect to ip T yields the above equation of motion ([5]) for the 
Majorana field. It should here be stressed again that the Majorana equation has 
two degrees of freedom less than the complex Dirac equation in its standard form 
and describes particle and antiparticlcs of opposite mean helicity. So using the 
real Dirac, i.e. the Majorana, equation implies and automatically ensures this 
reduction of the degrees of freedom (see the appendices Section[7]and Section[8]). 
Throughout the remainder of our paper we will stay with the real description. 
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3 Symmetries and Lorentz transformation in Ma- 
jorana representation 

Now let us briefly discuss the symmetries of the Majorana equation ([5]). To 
prepare this we recall that there is an important symmetry operator, namely 
the chiral matrix operator 7 s , which (see for example [2]) is defined as 7 s = 
i7°7 1 7 2 7 3 = ry 5 . Inserting the gamma matrices of equation ^ we obtain, 

7 5 = a x a y a z = f _° x * V (8) 

Like the other gamma matrices, 7 is purely imaginary. As a result one can not, 
by standard projection, obtain real right- or left-chiral components of the real 
eigen-spinor solutions of the Majorana equation (|5|) as given by (|42|) and (|43[) 
in the appendix Section [7J 

Obviously, the chirality operator 7 s by its definition (jSJ) commutes with a. 
But we require an operator like the r of Section [2] which anticommutes with the 
alphas. The product j3j 5 = 5 = (3a x a y a z has that desired property and thus 
corresponds to r. So we define: 



-7 
7 



(9) 



the square of which is -1 (like the square of r), and which anticommutes with 
/3 as well. 

Concerning the symmetries of the Majorana equation, we consider in partic- 
ular chirality conjugation C, parity V, and time reversal T operations. Generally 
speaking the Majorana equation is invariant under the symmetry operation O, 
if the spinor 

ip° = Oip (10) 

also fulfils that equation. Therefore, when applying the operation O from the 
left and its inverse O -1 from the right, whereby the unit operator is given by 
the decomposition OO^ 1 = 1, we obtain the result 

O— (IT 1 + Ola ■ ^O- 1 ~ mOfiO- 1 ) Ctyfx, t) = 0. (11) 
at ax / 

We conventionally define the time and space coordinate inversion operations T 
and P on a spinor ip by 

Tty(x,t)=tf(x,-t), (12) 
P^(x,i)=^(-x,i), (13) 

and also recall the complex conjugation operation C, which gives CiC -1 = 
— i, and but which here has no role to play as everything is real. With these 
preparations in mind, it is easy to see which operators provide the requested 
symmetry operations. We compose them in the Table [TJ To complete the 
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Table 1: Symmetry operations 



Operation Time reversal Parity Chirality conjugation 
Operator ST /3P 5 



operator algebra, it is important to note the following commutation relations. 
Of course, the coordinate reversal operators T and P commute with a, j3 and 5. 

Let us first consider in (jTTj) the time reversal, O = T = ST. Apparently, 
it affects the mass term via /3 as well as a, where the signs are reversed, and 
also changes the sign of the first term. Therefore, also ip 7 ' = Stp(x, —t) solves 
the Majorana equation ([5]). The parity operation O = V = (3P commutes with 
the mass term and does not affect the first term in ([5]), and it also leaves the 
momentum term invariant since a and x both change signs together. There- 
fore, also ip v = pip(—x,t) solves the Majorana equation. Finally, we consider 
chirality conjugation defined as O = C = S. It changes the signs of the alpha 
and beta terms in ([5]) with no net effect on the time-derivate term. Therefore, 
%jj c = Si/)(x,t) does not solve the Majorana equation, but in fact its chirality- 
conjugated version which is only derived in the appendix Section [5] In con- 
clusion, the symmetry operations of Table [T] work in a transparent way on the 
Majorana equation. 

To obtain the spin operator of the Majorana field, we now consider the 
Lorentz group generators, which in its four-component spinor representation 
are known [5] to be given by the commutator a^ v — g[7 ,i ,7 I/ ], which yields in 
Majorana representation with the gamma matrices as given in ([2]) the subse- 
quent matrix tensor: 



/ 





ia x 


ia y 






-ia x 











-ia y 









V 


-ia z 






o / 



(14) 



Note that these tensor elements are 4x4 matrices. Thus we obtained by defini- 
tion the spin operator S of the Majorana field, which unlike the spin operators 
in the Dirac or Weyl representation, is not diagonal and purely imaginary. It 
rCctds " 

E = i£ = i ((° ~o)'(l °)'(° "<?))■ <-) 

The spin operator does not commute with a but with f3. Therefore its 
cigcnfunctions are not eigenfunctions of the Majorana equation. Remember 
that a in turn anticommutes with /3. Note further that the spin operator is of 
course hermitian, and as a result we have S = — S, whereas a T = a. We 
also recall that the above spin matrices obey the same relations like the Pauli 
matrices, i.e. [S x , S y ] = 2iS z , where indices can be permuted in a cyclic way. 
So the Majorana spin operator can as usually be defined as Sm = ^S, which 
forms the SU(2) Lie algebra, also generating the familiar SU(2) Lie group of 
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the weak interactions. Finally, by using the a. and £ matrices, the Lorentz 
transformation operator A(i9, tp) turns out be real, and the related matrix can 
then be written as: 

S(A)=expQa.tf+i£.y), (16) 

where the vector angle i9 refers to genuine Lorentz transformations and ip to 
proper spatial rotations. Note that S(A) T j3S(A) = ft, a key feature which 
ensures Lorentz invariance of the Lagrangian, which is obvious for the mass 
term in particular in the Lagrange density (J7J. 



4 Gauge field theory of flavour dynamics 

We are now coming to the main theme of this paper, which is the question of 
Majorana (neutrino) gauge-field generated dynamics. Let us first consider an 
Abelian gauge field A^x) like in electrodynamics. Conventionally, this is intro- 
duced into the field equation by replacing, according to the minimal coupling 
principle, the time-space derivative by the covariant derivative 

D^ = d IM +ieA tl . (17) 

Concerning, the Majorana equation (|4]), this coupling causes a serious problem 
as the corresponding complex phase, like exp(ieA(cc)), where X(x) is a function 
of space-time and e the coupling constant (charge), in the wavefunction does 
not commute with the operator r. Fortunately, there is no such problem with 
(|5|), other than the solution spinor now has to be complex, which is not what 
we were aiming for to begin with. However, there is a remedy for getting rid of 
the unwanted imaginary unit, if we can find a non- Abelian symmetry group |22) 
the matrix representation of which is purely imaginary (and so electromagnetic 
interaction is excluded at the outset). The spin algebra as given by the Majo- 
rana spin operator according to (fT5|) provides just the desired representation, 
which yet is not identical with the fundamental unitary representation of SU(2). 
Therefore we will chose (with some coupling constant g) the connection involv- 
ing the un-normed isospin operator defined by the matrix three-vector S = iS 
as follows: 

D M = d ll + igS.A M = d ll -g8-A lt , (18) 

which is real and thus still permits real solutions to be obtained for the interact- 
ing Majorana field and its spinors. The real matrix isospin three-vector reads 
explicitly: 



S = 



/ 


f 











-1 


\ 


( 





-1 












( 





-1 







\ 










-1 







1 
























1 













1 





















-1 






1 















V 




1 











J 


V 








1 





J 






-1 








J 


) 



(19) 
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The spin matrix algebra leads for any real three- vectors A and B to the useful 
relation: 

(S • A)(S • B) = -A • B + S • (A x B). (20) 

As the Majorana matrix representation is four-dimensional, we have thus 
introduced in the Majorana field four new internal degrees of freedom related 
to the gauge field, which of course requires a physical interpretation. We return 
to this issue in the discussion section and here proceed formally. Note that S 
acts on a four-component spinor 



VV>4 J 



(21) 



Each single component relates to the Majorana equation (fTJ) or ([S]). For the 
free field the eigenfunction ip is given in the appendix Section [5] The extended 
Majorana Lagrange density, while including the gauge-field interaction, now has 
two parts which read: 



C M + A = i^ T (a^(d M - gS ■ A„) - 0m)¥. 



(22) 



We can identify in (|22[) a term associated with the particle flux density. This 
interaction term can conventionally be written as 



a = -r • a m , 



(23) 



where the real isospin current density generated by the SU(2)-like gauge sym- 
metry is given by 

= gi# T a"S* = ff * T a M S*. (24) 



The tensor of the gauge field strength is according to standard theory [221 H] 
given by the commutator of the connection (|18[) . i.e. we have 



F^u — — [Dfj, , Dy] — iS • F M „ , 



(25) 



where use has been made of (|20|) . Thus the antisymmetric field tensor (written 
as a normal three vector) is, again by help of (|20l) . derived in the concise form: 



c^A„ - d^Aft + 2gA^ x A„. 



(26) 



As is well known from non- Abelian gauge theory |221 [T] , the non-commuting 
nature of the vector components of the isospin S yields the nonlinear terms in 
the field tensor ([26]). The Lagrangian of the present SU(2)-like gauge field (see 
text books) theory is defined by 



C F = -l T r(F^FH = -1 (F^ • FH 



(27) 
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The term on the right of (j2"T)l comes from the properties of the isospin matrices 
in (Tl9l) . which also obey the relation: Tr(SiSj) = — 45jj-. This produces the 
normalization factor 16 in £p instead of 2 of the standard SU(2) gauge group 
in its fundamental representation. 

The two combined equations (|2"^j) and (|27p define the total Lagrangian, C = 
Cf + C\ + £mj and govern the dynamics of the fermionic Majorana field and 
the associated bosonic gauge fields. Variation of the Yang-Mills [15] action four- 
integral over C with respect to the relevant field variables finally yields the two 
real coupled partial-differential matrix equations: 

(t"0 m + m)¥(a;) = gfA^x) ■ Stffc), (28) 

d^F^ix) + 2g{A ll {x) x F»"(x)) = g^ T \x)6f 'Stf \x) . (29) 

Remember that 7^ = /3a M . Both equations (|28p and (|2"9")l have terms linear 
in the fields (free fields) and quadratic ones describing the field coupling. The 
gauge field equation (|29p in addition contains a term quadratic in this field, a 
property leading to nonlinear interaction of the field with itself. 

In the next section we will make some drastic approximations to the above 
complex system of coupled field equations, yet which provide an interesting 
simplified model which already seems to indicate the possibility of flavour mixing 
and oscillations at frequencies determined by the effective masses as induced by 
mean-gauge-field coupling. 

5 Flavour mixing and neutrino oscillations 

Let us consider a very much simplified version of the general gauge-field model 
equations ([28]) and (|29|) . Assume a fixed location and time dependence only, and 
further consider a mean-field approach to the gauge fields, thereby assuming a 
static field, with the vector A^{x) — (A, 0,0,0) being constant. Then j^A^ = 
/3A. Dimensionally, the coupling constant g times the field A must correspond 
to a mass, and thus we may write gA • S as a field-related mass matrix, which 
by help of (TT9J) takes the form: 



M 



( 


-My 


-M z 


-M x \ 


+My 





-M x 


+M Z 


+M Z 


+M X 





-My 


V +M X 


-AT. 


My 


J 



(30) 



which obeys M T = —M, where the superscript T indicates the transposed 
matrix, and which reflects the properties of the real Majorana field spin operator. 
Then the spinor time evolution is given by the mixing equation: 

/ S(^+M)+mW(t) = 0. (31) 
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In this simple mean-field model, there are four masses corresponding to the basic 
mass m and three others (inertia induced by mean-gauge- field coupling). The 
spinor solution involves oscillations or flavour mixing at the frequencies to be 
determined from the solution of the eigenvalue problem posed by equation ([3T]) . 
We can rewrite it as a standard second-order oscillation equation as follows. 
Multiplying the equations, i.e. their wavef unctions, by /? yields two coupled 
equations for and (3^, which can be inserted into each other to obtain a 
coupled oscillator (remember there are four flavour degrees of freedom) equation 
as follows: 

+ 2M^ - {Ml + M 2 + Ml) + m 2 ) = 0. (32) 

Without the field coupling we have independent harmonic oscillations with a 
phase mt, with it we obtain mixing, which is associated with frequency split- 
ting and damping or excitation in dependence upon the mean static gauge-field 
vector M = (Af x , M y , M z ). 

We may think that the mass term could be replaced with a diagonal mass 
matrix representing different intrinsic masses of the four flavours, i.e. we may 
define 

( mi \ 
m 2 
to 3 
\ m 4 J 

However, such a matrix fh does not commute with M, and the SU(2) flavour 
gauge symmetry would be broken by adopting it. Therefore, we will not consider 
this case any further for the sake of simplicity. But we consider now spatial 
variations, given again a constant background gauge field, and add the spatial 
derivative term to describe spatial propagation, and thus obtain: 

' 1 *(x,t) = ( J 5m + M)*(x,t), (34) 



(33) 



dt d: 



which differs from the free Majorana field ([5]) essentially by the gauge-field- 
induced mass mixing term. Considering again ^ and (3^!, and by eliminating 
one them, we obtain 



(S-£ +2 ^ + ^ 2+w2 ) vI/(x '^ - (35) 

Use has been made of the property (a ■ D) 2 = D 2 for any vector D, and that 
the alphas and beta anticommute. Since the equation is real, we can assume a 
harmonic wave superposition and make for the solution the ansatz: 

*(x, t) = £cos(£ , i - p ■ x) + 6 sin(£'t - p • x), (36) 

with four-component amplitude or flavour-polarisation spinors € and &. Inser- 
tion of this ansatz yields the coupled set of equations: 

(-E 2 + p 2 + m 2 - M 2 )€ - 2EM& = 0, (37) 
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(-E 2 + p 2 + m 2 - M 2 )6 + 2EM& = 0, 



(38) 



where M — \J — l/4Tr(Af 2 ) = VM 2 . For nontrivial solutions to exist the de- 
terminant of the system ([37| and (|38|) must vanish, which yields for the energy 
eigenvalues the four roots: 

E 1>2 = ±M + v / P 2 + m 2 , (39) 

E 3A = ±M - Vp 2 + m 2 . (40) 

For zero gauge field one retains the normal relativistic energies of free particles. 
Apparently, the constant gauge field yields energy splitting and flavour mixing, 
as the new eigenstates are obtained by mixing of the free-particle eigenstates. 
The mixed state is determined by the flavour spinors £ and ©, which we shall 
not calculate in detail. Owing to the application of the gauge field the four 
flavours acquired different energies in dependence on the field strength M, and 
thus the degeneracy of the original (just one mass m) energy spectrum has been 
lifted. 

Finally, we recall that for free Majorana particles of each flavour species 
the polarization spinors depend on energy and momentum, according to the 
eigenfunctions given in appendix Section [7J For the free Majorana field, when 
being written like in (|36|) in terms of sine and cosine functions, the corresponding 
four-component spinors C and S are determined by the polarization spinor a, i.e. 
C = l/\/2ReQ! and S = l/\/2lma for particles, and similarly for antiparticles. 
For them we just have to replace a by f3, whose dependence on energy is given 
by the coefficient e± as defined in the appendix Section [7] The corresponding 
flavour polarization spinors are more complicated, though, and have in principle 
to be determined by finding for each energy the eigenvectors of the full sixtecn- 
dimensional system given by equations (|37| and (|38|). This tedious calculation 
shall not be done here. 

6 Discussion and conclusion 

A gauge theory for massive neutral particles has been developed on the basis 
of the real four-component Majorana equation. The novel aspect is that by use 
of its spin operator a purely imaginary representation of the SU(2) algebra can 
be defined, which provides a covariant derivative or connection to the gauge 
field that is real. Such a minimal coupling can preserve the real nature of the 
Majorana equation. The associated isospin is four-dimensional, and thus via 
this procedure we introduce four intrinsic degrees of freedom to the Majorana 
field. What could be the nature of these degrees of freedom? 

The empirical fact that neutrinos oscillate and thereby change flavour has 
motivated us to make such a proposal. The main aim was to describe the 
mathematical feasibility of coupling the real Majorana field with a gauge field, 
a possibility that was not obvious and, to our best knowledge of the literature, 
not known before. Adding a fourth member to the flavour family, is here just a 
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matter of mathematical necessity and imposed by the choice of our isospin, but 
it remains physically speculative. However, flavour could in this way become a 
dynamic trait of the neutral Majorana particles, and thus lead to their linkage 
and a deeper understanding of neutrino mixing. 

We know the standard-model fermions come in three flavours, yet the present 
model would imply another fourth flavour, beyond the electron, muon and tau, 
and their neutrinos, and perhaps similarly for the quarks as well. The four 
flavours result here from the dimension of the representation of the gauge group, 
and as such are a cogent consequence if we accept that present description at 
all. Its dimension is that of space-time in the covariant Dirac and Majorana 
equations, and thus given just by the dimension of the gamma matrices. The 
particle quantum states in these equations are ordered according to past and 
future (antiparticles and particles) and right-helical or left-helical. The isospin, 
when being based on the spin operator that normally describes rotations in 
three-dimensional space, may correspond to the remaining four spatial orienta- 
tions, which is forward and backward and up and down. This may be considered 
a simple intuitive interpretation of the physical content of the 16-component 
spinor 'J, which arranges each flavour family into a quadruplett of real Majo- 
rana spinors for neutrinos, or perhaps complex Dirac spinors for charged leptons 
or quarks. So the sixteen components correspond to the sixteen (2 4 ) signed do- 
mains of space-time, and thus the number of flavours must be four but cannot 
be higher. 

The simple static mean-field approximation of the gauge field can produce 
Majorana-field flavour oscillations and mixing, can lift the mass degeneracy 
of the four Majorana neutrinos and cause a splitting of their mass spectrum, 
the size of which is determined by the gauge field strength. So mass is partly 
acquired here from gauge-field energy. This simple model just illustrates the 
physical potential of a flavour-mixing gauge theory of massive neutrinos. 

7 Appendix I: Eigenfunctions of the free Majo- 
rana equation 

Let us discuss briefly the eigenfunctions of the free Majorana field, which obeys 
equation ([3]) or (jj) and can be decomposed into its particle and antiparticle 
components, and according to the recent papers of Marsch [2U] be written 

V>(x,t) = ^p(x,t)+V>A(x,i). (41) 

These two contributions can be expressed, in terms of the complex four-component 
polarization spinors to be defined below, separately as follows: 

^p(x, t) = exp(— \Et + ip • x)a(p)&(p) + exp(iEt — ip • x)a*(p)a*(p), (42) 
V>a(x, t) = exp(-LBt + ip • x)6(p)/3(p) + exp(iEt - ip • x)6*(p)/3*(p). (43) 
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Apparently, i/'p.a = tpp aj an d thus the wave functions are real. The polarization 
spinors a and (3 are the two eigenfunctions of the helicity operator S • p (as 
defined in Section 13]) with eigenvalues +1 and — 1. 



(E • p)o(p) 

(s • pMp) 



-a(p). 



(44) 
(45) 



After quantization the complex Fourier amplitudes a*(p) and b*(jp) turn into 
creation operators of particles and antiparticles, which on average (in the sense 
of a quantum-field expectation value the Majorana field) have opposite helic- 
ities. For the sake of completeness we give here the full four-component po- 
larization spinors, which in terms of the angles of the momentum unit vector 
p = (sin 9 cos (f>, sin 9 sin </>, cos 9) and of the module p read as follows: 



a(p) 



(+£+(p) sin | e 
(+£+(p) cos | e" 



/ (+£+(p)cos§ e~20 + e _(p) s in| e+i^)(l - 

)(1- 
)(1 + 
)(1 + 



£-(p) cos | e 
e_(p) sin | e H 



/3(P) = . 



V (+e+(p)sin5 

' (— £_(p)sin' 
(+e_(p) cos : 



e+a? + e_(p) cos | e~ 

e ~i<P + e + (p) cos | e + 
c+i^ + £+(p) sin | e" 



(— £_(p)sin|e 2 — e + (p) cos | e H 



\ (+£_(p)cos§ e H 



£+ (p) sin | e" 



) 
) 

) / 

) 
) 

) / 



(46) 



(47) 



It was convenient to introduce the two real quantities 



£±{p) 



l E(p)±p 
2E(p) ' 



(48) 



1. The relativistic energy of a 



the squares of which add up to unity, e\ + e 
particle is E(p) — y m 2 + p 2 . A useful property of the epsilons is the obvious 
relation (E ± p)£j = me±. Using their properties, one can readily show that 
a* a = 1 and 13 ' (3 — 1, respectively, a) j3 = and (ft a — 0. If m — 0, then 
£ + = 1 and £_ = 0. It should be emphasized that the above spinors tAa and tpp 
always are a superposition of both helicity states, as it is required for a massive 
relativistic particle. The eigenvalue equation of the helicity operator (based 
on the spin operator, see Section [3]) in Fourier space reads for the complex 
conjugated polarization spinor: 

(E-p)a*(p) = -a*(p). (49) 

Similarly, for /3*(p), which has the opposite positive helicity, we obtain: 

(£.p)/3*(p) = +/r(p). (50) 

By multiplying out explicitly all possible scalar products involving the two com- 
plex conjugated polarization vectors a*(p) and j3*(p), one finds further that 
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&a* = {a T a)* = and ftp* = (/3 T /3)* = 0, but at/3* = (a T /3)* = \{e\ - e 2 _), 
and similarly fvot* = (/3 T a)* = i{s\ — e 2 _). Consequently, there are four inde- 
pendent real polarization vectors spanning the full Hilbert space of the Majorana 
equation. The eigenvalues ±E(p) are two-fold degenerate, and their subspaces 
are spanned by a(p) and a*(p), respectively $(p) and (3*(p). Both helicities are 
needed because of the finite mass of the relativistic Majorana particles. How- 
ever, the degeneracy concerning the replacement of the spin £ by its negative 
(spin reversal) has been lifted, because the real Majorana equation does not 
obey chirality conjugation (see again Section [3]). 



8 Appendix II: Two Majorana equations 

As was mentioned already in the main text, the two-component complex ([!]) 
and four-component real ([5]) Majorana equations are fully equivalent. However, 
as discussed by Marsch [THIHO], there is a second chirality-conjugated equation 
obtained by operating with r on (QJ, which yields 



d__ d_ 

dt 9x 



X(x,i) = -mrx(x,t), (51) 



with x = t~4>, which is the left-chiral counterpart to the right-chiral field <f>. This 
equation can be obtained by simply inverting the sign of the Pauli matrices 
in (j3J, including of course of r as it contains a y . The ± signs in front of the 
Pauli matrices reflect the reducibility of the Lorentz group, which can in fact 
be decomposed into its left- and right-chiral components. In the chiral or Weyl 
representation the charge conjugation operator is given by C = 7 V C = jt, which 
when being squared equals unity. So its eigenvalues are ±, and its eigenfunctions 
obey Ci/j± = ±"0±- The eigenfunction with positive eigenvalue obeys the Dirac 
equation and in Weyl representation is given by the spinor: 

*-(!>)-{**)■ 

Insertion of this spinor yields for the two two-component spinors <p and x the 
twin Majorana equations Q and (|51[) . which are coupled through the restricting 
condition x = T 4>i which according to (|52p guarantees chirality self-conjugation. 

Correspondingly, we obtain a second real four-component Majorana equation 
by simply taking the negatives of the matrices a, /3, 7 of equation (|3|) . This choice 
transposes the vector a of (|6|) to its negative and also gives a negative (3, so 
that we obtain the second real Majorana equation as 

d d \ ~ 

— - a ■ — 1 V(x,t) = -mPil>(x,t). (53) 

The solution is according to the Table Q] given by the chirally conjugated spinor 
ijj(x.,t) = 6ip(x, t), where t/j solves the Majorana equation (JSJ). Note that both 
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Majorana equations obey the parity and time-reversal symmetry, but individ- 
ually break by construction chirality conjugation which links them together. 
Concerning the spin operator (JTJ) , it will also be transposed for the second con- 
jugated Majorana field and turn into the negative of the spin of the primary 
field, i.e. S = — S. However, the connection to the gauge field after (fT5j) will not 
change, as once we have made our choice of the symmetry group their generators 
are given and do not undergo the discussed space-time symmetry operations. 

Returning to the twin two-component Majorana equations (U]) and (f5"Tj) . we 
recall that in the chiral Weyl representation of the Dirac equation the right- and 
left-chiral field are obtained by projection with the help of the chirality operator 
7 5 . Correspondingly, the Dirac field can be decomposed into its right- (index 
R) and left-chiral (index L) components, such that 

and (M) 

where we introduced the front factor for normalization, assuming that the func- 
tions and i/^ are individually normalized to unity. Above we used the 
symbols </>r = <f) and 0l = Xj with the chirality self-conjugation constraint that 
X = T(f) to be kept in mind. Therefore the mass term mixing the two chiral com- 
ponents, such that in the Lagrange density we have m-ifjip = m/2{(ftx + X^0), 
can formally be made diagonal by exploiting the above condition. The price 
to be payed for this decoupling is that the mass term becomes instead of triv- 
ial multiplication a nontrivial operator involving with r also the inconvenient 
complex-conjugation operation C. As the result we get the Lagrange density 

C = 2 (^R ( icr R 9 M ~ mr R) + ( icr L 5 M - mT L) <h) ■ ( 55 ) 

It can be shown (see the paper by Pal [23]) to be hermitian, while the two terms 
being their hermitian conjugates, which can be validated by using the relations 
= — t0l, respectively 0l = t0r, which are constitutive for the complex two- 
component Majorana field. Here we defined for the sake of formal symmetry the 
symbols: tr = r, tl = — t, with r = <7 y C, and <7p = (l,er) and cr^ = (L — <T )- 
However, with cr M = cr^ we may also rewrite (|55[) in terms of <p = 0R & s follows 

C = i Im (<jJ (icr^ - mr) 0) , (56) 

which emphasizes that there is only a single two-component complex Majorana 
field. 

The result of equation (|55|) can now be directly transferred to the real Majo- 
rana equation, following the same mathematics that lead to equation , respec- 
tively (|53p . By decomposing the two parts of (|55[) into their real and imaginary 
parts and writing the density in terms of four-component real spinors after [20] , 
the resulting alpha and beta matrices will of course change correspondingly, as 
we obtained qir = a and — —at, respectively, /3r = j3 and /3l = —ft- By 
definition, we obtain af^ = (1, ct) and = Consequently, af^ = <5l/j, 
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which shows how chirality and space-time are intrinsically connected. We can 
then write the combined Majorana Lagrange density as 

£ = \ (*r(or^ " /3r^)*R + ^lK^ - /3 L m)* L ) . (57) 

This equation shows a formal symmetry between the left-chiral and right-chiral 
components of the Majorana fields which at first sight seem to be independent. 
However, they are not, as we have = <^r, and vice versa ipR = —Sipi,. 
Also recall that {5, f3} = and {6, a.} — 0. Note that by its definition S T = 
—S and S 2 = — 1. So by inserting t/>l = fiipn into the above Lagrangian, the 
second term turns out to be exactly equal to the first. So there is only a 
single real Majorana field the Lagrangian of which was already given in ([7]). 
Therefore, considering the real Majorana field involving real four-component 
spinors, we believe, is advantageous over the complex Majorana field involving 
two-component complex spinors and requires to take care of the mathematical 
subtlety and complication of the operator r. 

9 Appendix III: Chiral decomposition of the spin 
operator 

Concerning the possible choice of the gauge symmetry group, traditionally 
SU(N) is used with the prominent SU(2) and SU(3) Lie groups employed in their 
fundamental representations for the weak and strong interactions. However, the 
real SU(2) representation chosen in (fT9|) is special in so far as it corresponds to 
the physical quantity angular momentum of the Majorana field. So there is no 
need to use left- or right-chiral Weyl fields, if the neutrino is assumed to be a 
massive Majorana particle, which by its very nature comes as a left- helical parti- 
cle and right-helical antiparticle. The Majorana spin operator Sm was adopted 
as isospin in the present work. It can be written in terms of the Pauli matrices 
in the form: 

s «=uu t)'(o- :).(: ~o°" )) ■ <*> 

We like to discuss then how to couple Dirac fermions (if also carrying the 
charge g) to the same gauge field as discussed in Section El The procedure 
how to do this is not at all obvious. Like in the weak interactions, perhaps a 
projection onto chiral eigenfunctions is needed for Dirac fermions. This is rather 
speculative and needs further investigations. As is well known the spin operator 
of the Dirac equation is given, both in Dirac and Weyl representation, by the 
expression: 

--K(o jl).(o- ;).(o- :,))■ m 

which is reducible and fully described by the fundamental SU(2) group as de- 
fined by the Pauli matrices. Here we want to mention that the spin can be 
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decomposed by means of the chiral operator 7 s and its related projection op- 
erators Pr,^ = 1/2(1 ± 7 5 ). They are idempotent, i.e. P^ L = Pr,l, provide 
a decomposition of unity: Pr + Pl = 1, and are orthogonal by construction: 
PrPl = 0. Furthermore, 7 s commutes with the spin operator, which is derived 
from the commutator of two gamma matrices after (|14[) . In fact it commutes 
with S in any representation. Therefore, we can write for the spin operator S, 
which obeys S x S = iS, formally 



and consequently obtain Sr : l x Sr.l = ISr^l, and [Sr, Sl] = 0. For the Majo- 
rana held the spin operator (TT5T) is purely imaginary and the projected compo- 
nents become complex, which is not helpful. In the chiral Weyl representation 
we simply get Sr : l = Pr,lc/2, which is irreducible and simplest. So the ques- 
tion comes up which the adequate representation is to be used for the SU(2) 
symmetry group. 

If we assume that any Dirac fermion is like a Majorana fermion endowed 
with four internal flavour degrees of freedom, then its interaction is mediated 
by the same gauge held A M considered for the Majorana fermion. We may take 
the isospin to be given consistently by the S of (TT9)) . and the coupling to the 
field be given by the same connection (fT5|) . However, we may prefer S = Sd to 
be used in the covariant derivative, i.e. favour the connection 



in the Dirac equation. Formally, such a gauge-held interaction model seems to 
resemble the weak interaction theory [16l |T7l [18] of the standard model, yet there 
are three major differences. First, when accepting the connection (|61[) for Dirac 
fermions, there is no restriction of the gauge-held-coupling to the left-chiral 
fermion field components only. Such a constraint need not be put anyway on a 
massive Majorona held, since it is chirally irreducible and by dehnition reduced 
to two degrees of freedom in comparison with a Dirac fermion. Secondly, both 
fermion species can be massive (with a single mass for all flavours), implying 
though that chiral symmetry is not obeyed. Thirdly, by construction of such a 
model, it will lead to lepton mixing and gauge-held mediated interactions among 
all fermions with different flavour, which are here assembled into unconstrained 
quadruplets, but not into doublets or singlets defined by chiral projection. 

Therefore, concerning possible choices of the isospin operator, we could make 
use of the individual spins as derived from the genuine but different gamma ma- 
trices in the Majorana and Dirac representations. Correspondingly, the covari- 
ant derivative would employ the respective spin operator of the held considered, 
which is then used as isospin providing the coupling to the common gauge held. 
This appears to be mathematically feasible, yet the physical implications remain 
unclear. To discuss these issues in depth is beyond the scope and intention of 
the present paper. 



S = (Pr + P l )S = Sr + S l , 



(60) 



Df> = dfj, + igS B ■ A 



(61) 
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